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Abstract. The quasi-isotropic inhomogeneous solution of the Einstein equations 
near a cosmological singularity in the form of a series expansion in the synchronous 
system of reference, first found by Lifshitz and Khalatnikov in 1960, is generalized 
to the case of a two-fluid cosmological model. This solution describes non-decreasing 
modes of adiabatic and isocurvature scalar perturbations and gravitational waves in 
the regime when deviations of a space-time metric from the homogeneous isotropic 
Friedmann-Robertson- Walker (FRW) background are large while locally measurable 
quantities like Riemann tensor components are still close to their FRW values. The 
general structure of the perturbation series is presented and the first coefficients of the 
series expansion for the metric tensor and the fluid energy densities and velocities are 
calculated explicitly. 



1. Introduction 

In the recent paper [I], we presented the generalization of the quasi-isotropic solution of 
the Einstein equations near a cosmological singularity found by Lifshitz and Khalatnikov 
[2] for the Universe filled by radiation with the equation of state p — | to the case of an 
arbitrary one-fluid cosmological model. Here, this solution is further generalized to the 
case of the Universe filled by two ideal barotropic fluids. 

As is well known, modern cosmology deals with many very different types of matter. 
In comparison with the old standard model of the hot Universe (the Big Bang), the 
situation has been dramatically changed, first, with the development of inflationary 
cosmological models which contain an inflaton effective scalar field or/and other exotic 
types of matter as an important ingredient |3J, and second, with the understanding 
that the main part of the non-relativistic matter in the present Universe is non- 
baryonic - cold dark matter (CDM). Furthermore, the appearance of brane and M- 
theory cosmological models jl] and the discovery of the cosmic acceleration j^j (see 
also jH] for a review) suggests that matter playing an essential role at different stages of 
cosmological evolution is multi-component generically, and these components may obey 
very different equations of state. Moreover, the very notion of the equation of state 
appears to be not fundamental; it has only a limited range of validity as compared to 
a more fundamental field-theoretical description. From this general point of view, the 
generalization of the quasi-isotropic solution to the case of two ideal barotropic fluids 
with constant but different p/e ratios seems to be a natural and important next logical 
step. The most straightforward candidates for these fluids are non-relativistic matter 
(CDM) with p = and radiation. One more popular candidate is ultra-stiff matter 
(p = e) which underlying physical model is a minimally coupled to gravity scalar field 
in the regime when its potential energy may be neglected as compared to its kinetic 
energy [Zj (in particular, it may be the inflaton field itself after the end of inflation). 

To explain the physical sense of the quasi-isotropic solution, let us remind that it 
represents the most generic spatially inhomogeneous generalization of the FRW space- 
time in which the space-time is locally FRW-like near the cosmological singularity t = 
(in particular, its Weyl tensor is much less than its Riemann tensor). On the other 



Quasi-isotropic solution near a cosmological singularity for a two-fluid model 3 

hand, generically it is very inhomogeneous globally and may have a very complicated 
spatial topology. As was shown in E] (see also |3]), such a solution contains 3 
arbitrary functions of space coordinates. From the FRW point of view, these 3 degrees of 
freedom represent the growing (non-decreasing in terms of metric perturbations) mode 
of adiabatic perturbations and the non-decreasing mode of gravitational waves (with two 
polarizations) in the case when deviations of a space-time metric from the FRW one 
are not small. So, the quasi-isotropic solution is not a generic solution of the Einstein 
equations with a barotropic fluid. Therefore, one should not expect this solution to arise 
in the course of generic gravitational collapse (in particular, inside a black hole event 
horizon). The generic solution near a space- like curvature singularity (for p < e) has a 
completely different structure consisting of the infinite sequence of anisotropic vacuum 
Kasner-like eras with space-dependent Kasner exponents |10| . 

For this reason, the quasi-isotropic solution had not attracted much interest for 
about twenty years. Its new life began after the development of successful inflationary 
models (i.e., with "graceful exit" from inflation) and the theory of generation of 
perturbations during inflation, because it had immediately become clear that generically 
(without fine tuning of initial conditions) scalar metric perturbations after the end of 
inflation remained small in a finite region of space which was much less than the whole 
causally connected space volume produced by inflation. It appears that the quasi- 
isotropic solution can be used for a global description of a part of space-time after 
inflation which belongs to "one post- inflationary universe". The latter is defined as a 
connected part of space-time where the hyper-surface t = t/(r) describing the moment 
when inflation ends is space-like and, therefore, can be made the surface of constant 
(zero) synchronous time by a coordinate transformation. This directly follows from 
the derivation of perturbations generated during inflation given in ^1] (see Eq. (17) 
of that paper) which is valid in case of large perturbations, too. Thus, when used in 
this context, the quasi-isotropic solution represents an intermediate asymptotic regime 
during expansion of the Universe after inflation. The synchronous time t appearing in 
it is the proper time since the end of inflation, and the region of validity of the solution 
is from t = up to a moment in future when spatial gradients become important. 
For sufficiently large scales, the latter moment may be rather late, even of the order or 
larger than the present age of the Universe. Note also the analogue of the quasi-isotropic 
solution before the end of inflation is given by the generic quasi-de Sitter solution found 
in [121 ■ Both solutions can be smoothly matched across the hypersurface of the end of 
inflation. 

Generalization of the quasi-isotropic solution to the case of multiple fluids provides 
a possibility to investigate the growing isocurvature mode of scalar perturbations and its 
effect on the growing mode of adiabatic perturbations in the regime when a space-time 
metric may not be globally represented as the FRW metric with small perturbations. 
Isocurvature perturbation of the energy density of a subdominant fluid need not be 
small, too. Note that the splitting of scalar perturbations into the adiabatic and the 
isocurvature mode is, to some extent, conventional (especially, in models with decay 
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of one forms of matter into others, e.g., in the so called "curvaton" model IE!)- Of 
course, this does not affect theoretical predictions for any observable quantity since 
contributions from all modes should be taken into account. We use the most natural 
splitting which is essentially the same as in the standard cosmological (A)CDM + 
radiation model, namely, that metric perturbations for the growing isocurvature mode 
are zero at t — 0. However, they are not zero for t > and grow with time. 

It is well known that the growing adiabatic mode remains constant in the 
synchronous reference system in the leading long-wave approximation for any one-fluid 
(and even any one-component matter) model both in the linear regime and with all 
back reaction effects taken into account. This immediately follows from the very form 
of the quasi-isotropic solution for one fluid - it corresponds to the factorization of spatial 
and temporal variables in its leading term. In a multi-component case, this property 
is known to be valid in the linear approximation, namely, that there always exists one 
solution of equations for perturbations (which we just call the growing adiabatic mode) 
which remains constant in the long- wave regime (see, e.g., [H| and for a more recent 
discussion J3|). However, some concern was expressed (e.g., in JE|) if this remains valid 
in the long-wave non-linear mult i- component regime. We will show that the factorization 
of variables in the first term occurs in the generalized quasi-isotropic solution, too, and 
this term remains the leading one either until spatial gradients become important, or 
until the energy density of a firstly subdominant fluid becomes comparable to that of 
the dominant one (that results in the change of the expansion law). 



2. Series expansion of the quasi-isotropic solution for a two-fluid case 

Let us consider a FRW model with two barotropic fluids satisfying the equations of 
state: 

Pi = he h I = 1,2 (1) 

where pi and E\ are the pressure and the energy density of the corresponding fluid. For 
usual fluids, < k\ < 1. However, similar to the case of one fluid, our solution remains 
valid for a larger range k\ > — ~ for both fluids. Let us take 

k 2 > fci. (2) 

As usually [HJ EI , we will work in the synchronous system of reference 

ds 2 = dt 2 — r ) a pdx OL dx^ . (3) 

Though this system of reference is not fixed uniquely (in other words, some of arbitrary 
functions entering into a solution written in this system are fictitious, or gauge artifacts), 
the choice of the initial cosmological singularity at the hypersurface t = in the quasi- 
isotropic solution (see Eq. ()20j) below) removes half of this ambiguity leaving only 3D 
spatial rotations in the leading term (with corresponding terms in higher orders) as the 
remaining gauge freedom. 
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To find the time dependence of the leading term of the quasi-isotropic solution, 
we write down the FRW energy equation in the vicinity of the cosmological singularity 
which reads 

+ ± = + A (4 , 

a 2 a 2 a Ql a" 2 

where A\ and A 2 are some constants, while the exponents a\ and a 2 are defined as 

ai = 3(l + fci), (5) 

a 2 = 3(1 + k 2 ) > a\. (6) 

We will look for a time dependence of the cosmological scale factor ait) (which 
determines the determinant of the spatial metric j a p) near the singularity in the 
following form: 

aW = E^ 7 "- (7) 

n=0 

Substituting this expansion into the Friedmann equation (@J) and comparing the smallest 
(i.e., most singular) powers in t, one can easily find that 

70 = — ■ (8) 

a 2 

Thus, the lowest exponent in the quasi-isotropic solution is defined by the stiffer fluid. 
The next term in the expansion ((7j) arises due to the presence of the second fluid in the 
Universe. Its time dependence is characterized by the exponent 

71 = 2 + 7o (l - ai) = 2 + (9) 

a 2 a 2 

The third term in the expansion (J7J) arises due to the presence of the curvature term. 
It has the same exponent as in the case of the one - fluid cosmological model pQ: 

72 = 2 - 7o = 2 - — . (10) 

a 2 

The next term of the expansion is again connected with the presence of the second fluid. 
The correspondent exponent is equal to 

2 4«i . , 

7 3 = 4 + 7o-27oai =4 + . (11) 

a 2 a 2 

Notice that three cases are possible depending on relation between a± and a 2 : 

1) 73 > 72, « 2 > 2«i - 2 , (12) 

e.g., if k 2 = 1, a 2 = 6 and k\ = 0, a\ = 3, i.e., the dominant fluid is ultra-stiff matter 
and the subdominant fluid is dust; 

2) 73 = 72, « 2 = 2ai - 2 , (13) 

e.g., in the case k 2 = 1, a 2 = 6 and k± = 1/3, a± = 4, so the first fluid is ultra-stiff 
matter and the second fluid is radiation, or in the dust-radiation case k 2 = 1/3, a 2 = 4 
and k\ = 0, a,\ — 3; 

3) 7s < 72, ol 2 < 2ai - 2 , (14) 
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(e.g., in the case k 2 = 1, a 2 = 6 and ki = 2/3, a.\ = 5). 

Now it is easy to understand that all terms which are present in the series expansion 
for the one-fluid quasi-isotropic solution remain in the solution for a two-fluid case 
because the mechanism of their appearance remains the same: interference between the 
dominant fluid and the curvature term. It is convenient to write a general formula for 
this set of terms in the following form: 

7n,o = 7o + n(2 - 270), n = 0, 1, . . . . (15) 

Another sequence of terms in the quasi-isotropic solution arises due to interference 
between the two fluids. It has the following set of exponents: 

7o,m = 7o + m(2 - 7 Qti), n = 1, . . . . (16) 

In addition, the full set of terms occurring in the quasi-isotropic solution includes terms 
arising from mixing between the two sequences (|T3j) and ()16|) . Thus, the expansion (|7jl 
can be represented as 

o(*)= E W 7 "' m (17) 



n,m=0 



where 



7n,m = 7o + ^(2 - 270) + m(2 - 7o«i). (18) 

It is easy to see that if the number 

2-2 7o = a 2 -2 = 3k 2 + l 
2 - 7o«i «2 - ax 3(/c 2 - 
is rational, a degeneracy between exponents 7„ im in Eq. (fTHj) is possible (an example 
of such degeneracy is given in Eq. (|13)l). However, it does not create any peculiarities 
in the structure of the quasi-isotropic solution. It is important that the expansion (|17|) 
with the exponents given in Eq. (|18|) contains all possible cross-terms. 

Correspondingly, the inhomogeneous quasi-isotropic solution for the spatial metric 
ja/3 can be represented as 

^ a/3 = ^ 7 (™' m )^ 2 70+n(2-27o)+m(2-7o«i)_ ^20) 
n,m=0 

For illustrative purposes, let us display the structure of this solution for some simple 
two- fluid models. For the model with radiation and dust, the series (|2~U|l can be rewritten 
as 

7^ = E7S^/2. (21) 

n=0 

For ultra-stiff matter and radiation, one has 

7^E7§i 2/3+W3 . (22) 

n=0 

In the case of ultra-stiff matter and dust, we get 

7<tf = 7S)* 2/3 + 7S* 5/3 + ffi + E 7$* 8/3+(n - 3)/3 , (23) 

n=3 
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and for stiff matter and the fluid with k% = 2/3, the solution reads 

7 Q /3 = E7^ 2/3+n/3 - (24) 

n=0 

Now, we will deduce explicit formulae for the first terms of the quasi-isotropic 
solution. The following notations will be used: 

Jaf) = a aP t 2 ^ + b a ^ 0+2 -^ + c a(3 t 2 + d aP t 2 ^ +4 ~ 2 ^ + ■■■. (25) 

As usually, the spatial metric a a p is arbitrary. For the time being we treat the metric 
b a/ 3 in the second term of the expansion (|2*5|) arising due to the presence of the second 
fluid as an arbitrary positively defined symmetric tensor. 
The inverse spatial metric reads 

_ c «/3^2-47o _ ^a/3^4-27o«i-27o _|_ yx^Sf3jA— 2700:1-270 {2Q) 

where a af3 is defined by the relation 

a Q % 7 = «f« (27) 
while the indices of all other matrices are lowered and raised by a a p and a a/3 , for example, 

b% = a a ^b l(3 . (28) 

Let us also write down the expressions for the extrinsic curvature, its contractions and 
its derivatives: 

Ka/3 = ^M = 2 l0 a af3 t 2 ^- 1 + (2 7o + 2 - 7oai)M 27o+1 " 7oQl 

+ 2c a(3 t + (2 7o + 4 - 2 lQai )d aP t 2 ^- 2 ^ , (29) 

«£ = 2705ft" 1 + (2 - joaMt 1 -^ + 2(1 - lo)^ 2 ^ 

+ 2(2 - 7ba 1 )d& 3 " 27ttai - (2 - 70^1) U^ 3-270 * 1 , (30) 

< = 670^ + (2 - 7o« 1 )6t 1 - 7oQl + 2(1 - joVj 1 ' 2 ^ 

+ 2(2 - 7o ai )^ 3 - 27oQl - (2 - 7o a 1 )6 a/3 6 a ' 3 t 3 - 27oai , (31) 



^ -2 7 o%t- 2 + (2 - 7oai )(l - 7o«i)^- 70 

+ 2(1 - 7o )(l - 2 7o )c^- 2 ^ + 2(2 - 7 o«!)(3 - 2 7o « 1 )^t 2 - 27oC 

- (2 - 7o ai) (3 - 27 a 1 )^ / ¥- 270Ql , (32) 

= -6 To t- 2 + (2 - 7 oax)(l - 7 oai)6^ 7oai 
+ 2(1 - 7o )(l - 2 7o )cr 2 ^ + 2(2 - 7o ai)(3 - 2 lo a 1 )dt 2 - 2 " / ° a ' 

- (2 - 7oai )(3 - 2 l0 ai)b aS b Sa t 2 ~ 2 ^ a \ (33) 

kIk% = 12 7o 2 t- 2 + 4 7o (2 - loai )br^ 

+ 8 7o (l - 7 o)ct- 270 + 8 7o (2 - 7o«i)rft 2 - 270Ql 

+ (2 - 7o ai)(2 - 70 ai - A l0 )b%t 2 - 2 ^ a \ (34) 
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Also, we need an explicit expression for the determinant of the spatial metric and its 
time derivative: 

7 = det lafi = t 670 det o(l + bt 2 -" 1001 + ct 2 ' 2 ^ 

+ dt 4 ~ 2 ^ ai + (1/2) (b 2 - b%)t 4 - 2 ^ oai ), (35) 

1 = fry,,*" 1 + (2 - TottO&t 1 - 7001 + 2(1 - To)^ 1 - 270 
7 

+ 2(2 - 7o a 1 )rft 3 - 270Ql - (2 - 7o a 1 )6^^ 3 - 27oai . (36) 
Now, using well-known expressions for the components of the Ricci tensor [T7j : 

(37) 
(38) 

;«2 (39) 

where P£ is the three-dimensional part of the Ricci tensor, and substituting the 
expressions f!29j) - (jSEJ) into Eqs. (|37|) - (p^5jl . one immediately obtains 

R° = 3 7o (l - 7o)t" 2 - ^(2 - 7o«i)(l - 7o «! + 2 l0 )bt~^ 



K = 






2 dt 


4 a ' 


R° a = 






R p = 


-Pf>- 








a 


2 <9t 



(1 - 7 o)ct- 270 - (2 - 7o«i)(3 - 27 q;i + 2 7o )dt 2 - 27oai 
1 

4< 



+ 1(2 - 7o«i)(4 - 3 7o «i + 4 7o )^6^ 2 - 27oai , (40) 



1 



< = i(2- 7 o«i)(6^-6 ;a )^ 7oai 



+ (1 - 7 o)«^ - c^t 1 - 270 + (2 - 7 o«i)(< /3 - rf ;a )t 3 " 270Ql 

+ ~(2 - 7o«i)((^); Q - (^^);/ 3 )^ 3 - 270Ql - (41) 

To write down an explicit expression for RP, we need an expression for pP taken up to 
linear terms in b a p. It reads 

p/3 — p/3y--27o _ yuv p/3 y.2-7 ai-27o 

+ \W;L «£ - C + W a f v )t 2 -^- 2 ^. (42) 

Here, P^ av denotes the three-dimensional Riemann-Christoffel tensor for the spatial 
metric a a p. Covariant derivatives are also taken with respect to this metric. Now 

ffP — _p/3£-27o _|_ fol*vp/3 ^.2-7001-270 

_ -(Oh v P — h P — h v P — \P v -I- h v & U2-70Q1-270 
2\ \va u ; a u ;au u a\v ' u a; v) 1 

+ 7o (l - 3 7o )5fr 2 - \{2 - 7o«i)((l - 7o «! + 3 7o )6f + 7oWf)t~ 7t,ai 
-(l-7o)((l-7o)c^ + 7oC^)r 270 
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- (2 - 7 o«i)((3 - 27 «! + 37o)d2 + 7o^)^" 27oai 
+ 1(2 - 7 o«i) ((3 - 270a! + 3 7o )« + 7o^£ 

-i(2- 7 o«i)^)t 2 - 2 ^. (43) 
The Einstein equations £1X6 , cLS usually, 

flj = 8nG(Ti - \5{T) (44) 
where the energy-momentum tensor for a two-fluid model has the form 

Tl = (ex + Pl )u u u lj - + (e 2 + p 2 )«2^ 2j - p 2 ^'. (45) 

Analyzing the Friedmann equation (J1J) again, one can see that series expansions for the 
energy densities (and, hence, for the pressures) have the following form: 

£ (n.™)^£-7oai+n(2-27o)+m(2- 7o ai) (45) 

n,m=0 

£ ^' m )^ x ^-2+n(2-27o)+m(2-7oai)^ (47) 

n,m=0 

It is easy to notice that the series expansion (|4"7|) contains all powers of t which appear in 
(|46|). However, the series (}4Tj) contains the term proportional to 1/t 2 which is absent in 
the series for e\. Besides, the series (jlTjl contains the terms £ ^ n ' ^t^ 2 + n ( 2 - 2 'yo) ^ n = 1,2, . . . 
which could appear in the series (jjpj) in the case of the degeneracy of the type given by 
Eq. (HU). 

Now let us write down components of the Einstein equations (|44jl in more detail: 
R° = 4vrG((l + 3Jfei)£! - 2(1 + fci)£i«i a < 

+ (l + 3fc 2 )e 2 -2(l + A; 2 )e 2 ui a i£) J (48) 
i2° = 8ttG((1 + kjewxaul + (1 + A: 2 )£ 2 m 2q m°), (49) 

Ri = 4ttG(2(1 + fcOeiUteW? - (1 - fc^e^ 

+ 2(1 + A; 2 ) £2M2aM f - (1 - fc 2 )e 2 0- (50) 



Comparing the structure of the left-hand side of Eq. (J49J) with its right-hand side, 
one can see that the series structure for spatial components of fluid four-velocities u a 
has the following form: 

«2a= E ut m) t 1+n{2 ~" /0)+m(2 ^° ai) (51) 
n,m=0 



where w 2 a°^ — 0, and 



Ula = J2 M ^tl+«(2-70)+ m (2-70ai). ( 5 2) 
n,m=0 

Correspondingly, 

M « _ ^ M a(n,m)^l-27o+n(2-7o)+m(2-7 ai) (53) 
n,m=0 
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where tio' 0,0 ' = 0, and 



u ~ _ u ^ n ' m ^ 1 - 2 70+n(2-7o)+m(2- 7o a 1 )_ 



■" — ^ t , a ( n ' m )fl-270+"(2-7o)+m(2-7oai) 
n,m=0 

Also, we need the following expressions: 

U 20 = U° 2 = 1+ Yl ^-)i 2 - 2 T0+n(2-70)+m(2- 70 a 1 ) (55) 
n,m=0 

where i4o°' > = 0, and 

UlQ = u = 1+ M ^ m )t2-270+n(2- 70 ) + m(2-70 Ql ). (55) 

n,m=0 

Now, taking the leading term proportional to t~ 2 in the expansion for R® and comparing 
it with the corresponding terms in the series for the energy-momentum tensor, one can 
find from Eq. f48|) that the leading term in the series for the density energy (|4"?|) is given 
by the formula 

(0,0) _ 1 ( z 7) 

Thus, the firstly dominant (stiffer) component becomes homogeneous at t — > 0. 

The next step consists in comparison of the terms proportional to t~ 2 ~ t ° ai in the 
equations ()48j) and (|50|) . As a result, two consistency relations for the two unknown 
variables e^' 1 ^ and arise: 

- (2 - 7o«i)(l - 7o«i + 270)6 

= 8ttG((1 + 3h)ef' 0) + (1 + 3A; 2 )4°' 0) , (58) 
(2 - 7 o«i)((l - 7oai + 3 7o )&f + 

= 8ttG((1 - fei)e?' 0) ^ + (1 - fc 2 )4°' 1) <5f). (59) 
It follows from Eq. ()59|1 that 

bi = \b€- (60) 

This means that we cannot choose the metric b a @ freely; it should be proportional to 
the metric a Q( g, i.e., 

b a/ 3 = l^ba a p. (61) 

Thus, only one new arbitrary function of spatial coordinates b(r) appears in the second 
term of the series in the right-hand side of Eq. (J25|) . We will immediately see that 
it describes the amplitude of isocurvature perturbations. Really, substituting the 
expression (JHT|) into Eq. (jBUjl . one get 

(2 - 7o«i)(l - 7 «i + 670)6 

= 247rG((l - kx)ef ,0) + (1 - A; 2 )4 0ll) ) (62) 
and then, resolving Eqs. (|55jl and (p32"|) with respect to ef^ and £ 2 , we obtain 

(0j0 ) _ (3fc 2 - 2fci + 1)6 
£l ~ 12nG(k 2 + 1) 2 ' [b6) 
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*" - -u^ohrr <64) 

Therefore, the firstly subdominant component 1 is generically essentially inhomogeneous 
near the singularity 

Now, we should compare terms proportional to t~ 2l ° in the left- and right-hand sides 
of Eqs. (jlHj) and (J57Jj) . We consider a generic situation when the degeneracy described 
above (see Eq. ()13|)) is absent. Then we have 

-(l-7o)c = 47rG(l + 3A; 2 )4 1 ' 0) , (65) 
^ + (1-7o)((1 + 7qK + 7o^) 

= 4vrG(l - fc 2 )£?' 0) 5f. (66) 

Resolving Eqs. (JHSJ), (|6l)J) with respect to variables and e^^, one obtains the 
following expressions which coincide with known expressions for the one-fluid quasi- 
isotropic solution pQ: 

* 9(fc 2 + l) 2 ( pi3 P(3fc 2 2 -6fc 2 -5)(% \ 

a " (3fc 2 + l)(3A: 2 + 5) V " 4(9fc 2 + 5) 1 J 

r d,o) _ 3P(fc 2 + 1) 

£2 " 16ttG(9A; 2 + 5)- m 

The next step consists in the consideration of terms proportional to t 2 ~ 2loOL1 . Using 
Eq. (|6T)|) . one can obtain from Eqs. PH|) and (jBTI|) that 

(0 wo o w i (2-7oQi)(4-37oa 1 +4 7o )6 2 
- (2 - 7o«i) (3 - 2 7o ai + 2 7o )(i H — 

= 47rG(4°' 2) (l + 3A; 2 ) + ef'^l + 3A*)), (69) 

(2 - 7 o«i)((3 - 2 7 o«! + 3 7o )^ + 7 oO 

(2-7 ai)(12-77oai + 127o)6 2 ^ 
36 

= 4vrG^(4°' 1) (l - fei) +4°' 2) (1 - k 2 )). (70) 
It is easy to see from Eq. (fTOj) that 

€ = \d%. (71) 
Substituting Eq. (JTT|) into Eq. (JTUJ) , one get 

(2 - 7q«i)(3 - 270^1 + 6j )d5P 
3 

(2 - 7o a 1 )(12- 770^ + 1270)^ 
36 

47rG5f(£S ' 1) (l - h) +4°' 2) (1 - fe)). (72) 



+ 
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Now, using the definitions of the coefficients cti, a 2 and 70, it is convenient to rewrite 
Eqs. (JoTJj) and (|T2*|) in the following form: 

e { 2 °' 2 \l + 3k 2 )+ef' 1 \l + 3k 1 ) 

[k 2 - h) (9k 2 - 12k, + l)d (h-k^Gh-dh + l) 



6nG(k 2 + l) 2 3QnG(k 2 + l) 2 

{k 2 -k l ){3k 2 -Ak l + 3)d 

QnG(k 2 + l) 2 367rG(A; 2 + 1) 2 ' 1 ' 

The system of equations (f73^) and (f74"j) contains three variables d, e^' 2 ^ and ef' . 
To solve it we use the energy conservation law 

^0 ;i = 0. (75) 

In the case of non-interacting perfect fluids, Eq. (|75]l gives two separate equations in 
our approximation: 

£2 + ~(£2+P2) = 0; (76) 
2 7 

£i + ™(ei+Pi) = 0. (77) 
27 



Using the series f)4fij) . (|47jl and (|HEJ), we get the following formulae connecting 

(0,2) _(0 
2 5 fc l 



£, ',£ ' , and b: 



-(0,2) _ j (3fc 2 + 5)6 2 

2 ~ 12ttG(1 + fc 2 ) 1447rG(l + Jfe 2 ) ' 1 J 

(0,1) _ (fc 1 + l)(3fc 2 -2A; 1 + l)6 2 
£l 24ttG(1 + A; 2 ) 2 ' 1 } 

Substituting expressions for the energy densities from Eqs. (|7Sj) and (J79|) into Eq. (J73J), 
one get the following expression for d: 

, (-3k 2 + 6k 2 - 6A;iA; 2 - 2A; 2 - 2k x - l)b 2 . , 

?■ 1 ; , , 2 — '—. 80 

12(5Jfe 2 -4A;i + l) v ; 

It is easy to check that this expression satisfies Eq. too. Substituting (J80|) into Eq. 
(1781 . we have 



(3^ - kj - hk 2 + 5k 2 - 3h + l)b 2 
1 " 24ttG(1 + A; 2 )(5A; 2 -4A; 1 + 1) 



-(0,2) 



'I) 



Now, we would like to write down explicit formulae for first coefficients of the 
velocity series (JoTj) . (J52J). Using these series and the series for the energy densities (JH)j), 
(|47|) . one can get the following relation for terms proportional to t 1_ Toai j n ^he Einstein 
equation (jJOJ): 

i(2 - 7 o«i) (62 ;/3 - 6. a ) = 8ttG((1 + fc 2 )£| ' 0) 4 / ) + (1 + *i)e£ 0) «£ 0) ).(82) 
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Eq. (}82|) contains two unknown quantities u^Jf 1 and uf^- Thus, we need additional 
relations to resolve it. We will use remaining components of the energy- momentum 
conservation law: 

T a ;i = 0. (83) 
In the synchronous system of reference, this equation can be rewritten as 

T° a , + 7? # + \K T l = 0- (84) 
Using the definition of energy-momentum tensor components for perfect fluids and 
the explicit expression for the trace of the extrinsic curvature ([31)1. we have the following 
couple of relations in the needed order of perturbation theory: 



he^l = (1 + k^er^) + ^(1 + k 2 )er ] u^\ (85) 

hef'l = (1 + h)§- t (sruT) + ^(1 + h)sruT. (86) 

Substituting the expressions (JH7|) . (jMJ) and (jHHj) for the coefficients £j ° , 4 and 4°'°' ) 
into Eqs. and (|HB|l. one finds 

(0,1) _ fc 2 (fc 2 + l)6;q , , 

-"2(^-2^ + 1)' (87) 

(o,o) _ fci(fc 2 + 1)6 ;Q 

+ l){h 2 -2k 1 + 1)6' 

The direct check shows that the velocity coefficients (j87jl and (}88j) satisfy Eq. (|82*|). 

The next order term of the quasi-isotropic series solution (proportional to t 1-2 T° in 
the Einstein equation (|49p) immediately gives 

(1 - 7o )(cJ ;/3 - c. Q ) = 8vrG(l + k 2 )e$ fi) v& 0) . (89) 

Notice that in the absence of degeneracy between exponents of different terms in the 
perturbation series, this equation contains the contribution from the firstly dominant 
fluid 2 only. Thus, using the explicit expression for the coefficients c£ (see Eq. (IHTjl and 
the Bianchi identity for the curvature tensor 

Pa = \P,a , (90) 

one has 

( o,o) _ 27k 2 (k 2 + lfR a 

2a 8(3jfe 2 + 5)(9A; 2 + 5)" { ] 

To find next terms of the quasi-isotropic solution for the velocities, i.e., corrections 
proportional to £( 3_7oai ), we shall use the identity (|84|) again. This leads to the following 
relations : 

k 2 efl = (3 - 2 7oai )(4°' 0) 4° Q 2) + e^u^) 

+ ^ 7 o(4°' 0) ^ 2) + 4°' 1) ^ 1) ) 

+ i(2- 7oai )fer ) 4 « 1) , (92) 
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he?S = (3 - 2 7o a 1 )(4°' C M Q 1) + s^uT) 

+^ 7o (,r ) ^ i) +4 o ' i) 4r ) ) 

+ i(2- 7oai )fer ) ^ 0) - (93) 

From Eqs. (j^j) . (J9l?j) one can find the formulae for unknown quantities and 
which read 

M ^ 2) = 7^ — au^u m x {h{k2 + 

(3fc 2 - 4fti + 1)4 

- (3A; 2 - 4fe + l)^ 1 ^ - (k 2 - h)be?>%^), (94) 

= i^i — au m x ^ h + 1)£ i°5 

(3fc 2 - 4fci + ljei 

- (3k 2 - Ah + l)^ 1 ^ - (k 2 - h)be?>%?f). (95) 

Substituting the explicit expressions for ef'°\ef' 1 \ef' 2 \ef' \ef' 1 \uf^ 1 \uf ( f^ from 
Eqs. flEZI), (El, (EH), dZSD, (E3), (JHi into Eqs. JMD and flHSJ), we get 

(0,2) k 2 (k 2 + l)bb ]a 

U2a ~ 4(k 2 - 2h + l)(5fc 2 - 4/n + l)(3/c 2 - 4A;i + 1) 
x (-9^ + 18A;i^ - 14A; 2 £; 2 + 4Jfe? + 3fc| 

- 6/fc — lfc 2 + - l 2 + 5fc 2 — 4A; - 1 + 1), (96) 
(o,i) = fcg(fc 2 + l)(l-fc 1 )(fc 2 -2fc 1 -l)& ;a 

UlQ 2(fc! + l)(3fe - 4Jb! + l)(Ar 2 - 2fe! + 1) ' 1 ' 

In conclusion, let us write down the expressions for the first terms of the quasi- 
isotropic series solution for the metric, energy densities of two fluids and their 3- velocities 
calculated above: 

4 2(3fc 2 -3 fc l+ 2 ) 

7a/3 = a a pt^^ + ba a(3 t 

9(/c 2 + l) 2 (~ ( (3£; 2 -6£; 2 -5)Pa a/3 \ 2 



12(5A: 2 -4fci + 1) 

1 6 3(fci+l) 
£9 = £ fc 2+l 

6vrG'(A: 2 + l) 2 t 2 127rG(A; 2 + 1) 
3(k 2 + l)P 



+ 4(9A; 2 + 5* !/ 

4- Ok -U1 U 2 n „ 

t 3 ( fc 2 +1 ) + • • • , (98) 



(3fc 2 + l)(3fc 2 + 5) V " p 4(9A; 2 + 5) 

(3fcf - 6A; 2 + 6fcifc 2 + 2£; 2 + 2fci + l)& 2 a a/3 2(6*2-6^+2) 



j ; — _ £ 1 3 < fc 2+i) 

16vrG(9A; 2 + 5) 

(3A; 2 - kj - hh + hk 2 - 3fci + ljg gfe-^pl) 

24vrG(A; 2 + l) (5^-4^ + 1) ' 1 j 

(3fc 2 - 2*! + l)b 2 (fcl +i) 

£1 = - — — — —t fc 2+ 1 

1 12ttG(A; 2 + l) 2 

(jfc 1 + i)(3fc 2 -2fc 1 + i)6 2 ga^za 

247rG(A; 2 + 1) 2 ' 1 J 
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k 2 (k 2 + l)h a ■ 3 * 2 -a*i+i 

27* a (k + 1)*P ;Q 



+ 



(3A; 2 + 5)(9& 2 + 5) 

fc 2 (/c 2 + l)66. a 



4(fc 2 - 2k! + l)(5fc 2 - Ah + l)(3fc 2 - 4A;! + 1) 
x (-9Jfcf + 18A;i^ - 14Jfe?Jfc 2 + 4fc? + 3^ 

- 6A;iA; 2 + 2k\ + 5A; 2 - Ak x + l)t k * +1 + • • • , (101) 

= h(k 2 + l)b, Q 
Ula (hi + l)(k 2 - 2ki + 1)6 

, fc;(fc 2 + i)(i-fc 1 )(fc 2 -2fc 1 -i)6 ;a a^±i fl02l 

2(fci + l)(3fc 2 - 4fci + l)(fc 2 -2fci + l) ' 1 J 

Note that the velocity flow is potential. Actually, it can be shown that this 
important property remains in all higher orders of the quasi-isotropic solution. In the 
degenerate case (see Eq. (fT3|)) when the relation between coefficients k\ and k 2 is 

h = (103) 
o 



the terms with superscripts (1,0) and (0,2) in the expressions (JHEJ), (j99|l and (jlOljl . 
i.e., the terms proportional to P and b 2 respectively belong to the same order in the 
quasi-isotropic series expansion. Then the coefficients at the corresponding powers of 
time t should be simply added due to linearity of all iterative equations used above. 



3. Conclusions and discussion 



The equations (|98 p - ()102jl represent the main result of the present paper - the first terms 
in the generalized quasi-isotropic series solution for the case of two barotropic fluids. 
How many arbitrary physical (i.e., gauge invariant) functions of spatial coordinates 
does it contain? The answer is four: 3 functions are contained in a a/ g (no conditions 
are imposed on this symmetric tensor, and the remaining freedom of coordinate 
transformations at t = consists of 3 spatial rotations ) and the fourth function is b. The 
former 3 physical functions describe the growing mode of scalar adiabatic perturbations 
and the non- decreasing mode of gravitational waves (with 2 polarization states) in the 
regime when metric perturbations are not small, while the latter function describes the 
growing mode of scalar isocurvature fluctuations. 

It is seen that the first term of the generalized quasi-isotropic solution remains 
factorized with respect to t and spatial coordinates. This represents the generalization 
of the known result, that small relative metric perturbations in the synchronous system 
of reference corresponding to the growing adiabatic mode remain constant in the long- 
wave limit even in a multi-fluid case, to the case when metric perturbations are not 
small. For the other, isocurvature mode, metric perturbations are zero at t = 0, but 
grows with time. However, their back reaction appears in higher order terms of the 
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generalized quasi-isotropic solution and does not affect the leading term of the adiabatic 
mode (as far as the series expansion has sense at all). 

However, as compared to the one-fluid case, the following new feature arises. In 
that case, the quasi-isotropic series solution loses sense in the course of future expansion 
when spatial gradients of all quantities become comparable to temporal ones. The 
series for the generalized quasi-isotropic solution may diverge for another reason even if 
spatial derivatives are still small, namely, when higher-order terms containing b become 
of the order of the leading term. This means that the firstly subdominant fluid becomes 
dominant and vice versa. Then the law of expansion should change, too (as occurs, e.g., 
in the standard CDM+radiation FRW model after the moment of matter-radiation 
equality). This shows that all terms in the generalized quasi-isotropic series containing 
the function 6(r) but not its derivatives can be summed. We leave this question for 
future work. 
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